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Abstract 

We consider power-behaved contributions to hard processes in QCD arising 
p • from non-perturbative effects at low scales which can be described by intro- 

p • ducing the notion of an infrared-finite effective coupling. Our method is based 

D . on a dispersive treatment which embodies running coupling effects in all or- 

^ ' ders. The resulting power behaviour is consistent with expectations based 

on the operator product expansion, but our approach is more widely applica- 
^ ■ ble. The dispersively-generated power contributions to different observables 

■ are given by (log-)moment integrals of a universal low-scale effective coupling, 

with process-dependent powers and coefficients. We analyse a wide variety of 
quark-dominated processes and observables, and show how the power contri- 
butions are specified in lowest order by the behaviour of one-loop Feynman 
diagrams containing a gluon of small virtual mass. We discuss both collinear 
safe observables (such as the e~^e~ total cross section and r hadronic width, 
DIS sum rules, e~^e~ event shape variables and the Drell-Yan K-iactor) and 
collinear divergent quantities (such as DIS structure functions, e~^e~ fragmen- 
tation functions and the Drell-Yan cross section). 
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1 Introduction 



Power-behaved contributions to hard colhsion observables are by now widely recognized 
both as a serious difficulty in improving the precision of tests of perturbative QCD and as a 
way to explore non-perturbative effects. Such contributions are manifest as power-varying 
discrepancies between fixed-order perturbative predictions and experiment. For some 
quantities the leading power contributions are of order A/Q where Q is the relevant 

energy scale and A ~ 300 MeV is the characteristic scale of QCD. In this case, even at 
scales Q ~ Mz ^ A, power-behaved terms can be comparable to O (a^) perturbative 
contributions. In other processes such as deep inelastic lepton scattering (DIS) and 

Drell-Yan lepton pair production, contributions of order /Q"^ are important because Q 
is not so large. At even lower scales, as in the hadronic decay of the r lepton, a stronger 
suppression of power contributions is needed in order for a perturbative analysis to be 
applicable [P-[T3[|. 

Such contributions are usually regarded as "power-suppressed corrections" to the per- 
turbative prediction, in which context they might appear negligible in comparison with 
perturbative terms of any order. From this viewpoint the former could not be considered 
before summing up the entire perturbative series, which in itself is ambiguous. In practice, 
however, power contributions can easily be distinguished since they are rapidly varying 
in contrast to the slowly (logarithmically) varying fixed-order perturbative contributions. 
Thus we can rather regard them as power enhanced, in the sense that they grow more 
rapidly as the hard process scale Q decreases. In view of their phenomenological rele- 
vance, it is important to study any possible source of power-behaved terms which might 
be a general feature of non-perturbative QCD. 

So far, there is no rigorous systematic method for calculating power contributions in 
QCD. In some cases the operator product expansion (OPE) provides a framework for the 
classification of possible power-behaved terms, but it gives little indication of their coeffi- 



cients ll^Jlql . General studies of the divergence of the perturbation series in high orders 
have led to the concept of renormalons |16| , which are singularities in the complex plane of 
the Borel variable conjugate to the running coupling. In particular, infrared renormalons 
lie on the integration contour of the Borel variable. In regularizing these singularities one 
generates power-behaved contributions whose coefficients are ambiguous. Formal mathe- 
matical manipulations alone cannot resolve the infrared renormalon problem, which is of 
a physical nature, namely the notion of the running coupling in the small fc^ regime. 

In this context, it is natural to consider the possibility that a finite running coupling 
might be definable at low scales, at least as an effective measure of the strength of inter- 
action at large distances. One may hope in this way to take into account the main effects 
of confinement at a sufficiently inclusive level. This idea can only be meaningful if the ef- 
fective coupling so introduced is found to be universal in some approximation. Its specific 
form, aefr(fc^), at small k"^ then determines the power contributions. Thus the experi- 
mental analysis of these contributions will constrain the form of aeff(/c^) at small scales. 
The theorists' task is to analyse the power contributions which are generated by the non- 
perturbative regime of the running coupling. To this end, one should study perturbative 
Feynman diagrams for hard process observables, in which the scale of the running coupling 
is not restricted to the hard region, but also runs into the non-perturbative domain. 

A technique for studying power contributions within the language of perturbative 
QCD is to introduce a small gluon mass p|J^, p!^ , p!^ , which may be used as a "trigger" 
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for long-distance contributions. In the present paper, we attempt to put this approach 
on a somewhat clearer theoretical basis. We work with a gluon field of zero mass, and 
introduce a dispersive representation of the running coupling ^^(fc^), which is assumed to 
remain finite all the way down to k"^ = 0. In this way ^^(/c^) is expressed in terms of an 
effective coupling aeff(/i^) depending on the dispersion variable /i. In this representation /x 
plays a role similar to a gluon mass as far as the phase space and Feynman denominators 
are concerned. The region of small values of fi is responsible for power contributions. 
Since the running coupling is by assumption universal, these contributions can be related 
for different processes. We study the corrections to a wide range of quark-dominated 
processes due to a single off-shell gluon with time-like or space-like virtuality. The results 
should provide the means to test experimentally the concept of a universal infrared-finite 
coupling as an effective measure of strong interactions at large distances, as reflected in 
sufficiently inclusive hard process observables. Our discussion is mainly in the context 
of Abelian gauge theory, although we argue that its application to the non-Abelian case 
may be justifled for leading power contributions. 

In characterizing the effective coupling at low momentum scale we assume, in agree- 
ment with the standard ITEP-OPE approach [|l5l , that non-perturbative physics does not 



affect the high-momentum region, that is, the propagation of quarks and gluons at short 
distances (soft conflnement). Within our method, this implies that the non-perturbative 
"effective coupling modiflcation" , (5aefr(At^), which is responsible for power corrections in 
hard distributions, is essentially restricted to the small /z^ region and must not directly 
introduce power corrections into as itself (notice again that fi^ is a dispersive variable 
and we are working in Minkowski space). To satisfy this condition, Saesif^"^) must have 
vanishing moments with respect to /i^, at least for the first few (integer) powers of /i^. 
If this were not the case, the predictions of the operator product expansion would be 
invalidated. 

We start in Section 2 by discussing the extent to which a running coupling might be 
definable beyond perturbation theory. We introduce a dispersive representation of as in 
terms of a spectral density function Ps(/i^) and then in terms of the effective coupling 
«efr(/^^)- We recall how this representation leads to a natural choice of the scale of as 
in deep inelastic scattering, and we define the corresponding 'physical scheme' which we 
adopt for the QCD coupling in the perturbative region. Next we define the effective 
coupling aeff(/i^), related to the spectral density function, which will be a central quantity 
in our treatment of power contributions. 

In Section 3 we present the method by which we compute leading power contributions 
from integrals involving the "effective coupling modification" 6aes- We show that in 
quark-dominated processes, ie those involving no gluons at the Born level, the leading 
effect of the running coupling is determined by the behaviour of aeflf(/i^) and a process- 
dependent characteristic function T . This is a function of the ratio e = ji^ /Q"^ , together 
with any relevant dimensionless ratios of hard scales {x}, given by the sum of all one- loop 
graphs computed with a non-zero gluon mass /i. The power contributions that we seek 
to estimate are associated with the non-analytic terms in the small-e expansion of this 
function. Each such term implies a power correction proportional to an integral of 5aeff(/^^) 
times a corresponding non-analytic weight function. Thus the task becomes to compute 
jF(e, {x}) for various processes and to see whether the predicted power contributions are 
consistent with experiment for some choice of the behaviour of the "effective coupling 
modification" 5aes at low scales. 



In Section 4 we present results on the characteristic functions for a variety of processes 
and discuss their behaviour at small e and the associated power contributions. In many 
cases we find strong enhancement factors, either logarithms of or singular functions 
of the auxiliary variables {x}, which lead us to hope that these contributions may be the 
dominant non-perturbative effects, at least in some important regions of phase space. In 
the case of A/Q contributions, no other possible sources of such terms are known, and 
therefore our results should provide useful semi- quantitative estimates of their relative 
magnitudes in different observables. 

We do not undertake any detailed phenomenological analyses in the present paper, 
but confine ourselves to providing the theoretical results on which future analyses could 
be based. We end with a summary of results and a discussion of the limitations of our 
approach and its possible further extensions and applications. 



2 Running coupling 

In hard processes initiated by a quark, the inclusive quantities at a given order in pertur- 
bation theory receive contributions from an extra gluon both from virtual corrections to 
the amplitude at the previous order (fc^ < 0) and from a new production channel {k"^ > 0). 
We discuss here how one can associate a running coupling with this gluon. For large vir- 
tuality \k'^\ it is clear how to identify the running coupling and the proper scale for its 
argument using perturbative methods. For small momenta, k"^ —>■ 0, the very language of 
quarks and gluon becomes scarcely applicable, and the notion of the running QCD cou- 
pling becomes more elusive. This situation is opposite to the case of the Abelian theory, 
in which the physical fine structure coupling is defined at vanishing photon momentum, 
a{0), and appears, for example, in the Thomson cross section. 

We explore here the hypothesis that in a non- Abelian theory the notion of the running 
coupling can also be extended to the region of vanishing momenta, at least in some effective 
sense. In doing so we seek to use the language and methods of perturbative QCD to probe 
some non-perturbative phenomena, in particular the leading power contributions to hard 
cross sections. In order to formulate this hypothesis more precisely, we assume that in 
QCD the running coupling can be represented by a dispersion relation which is inspired 
by the Abelian theory. Therefore we first recall the case of QED and then describe the 
extension to QCD. 



2.1 Space-like gluon 

Abelian case. Consider first the exchange of a photon with negative virtuality — fc^ 
(hereafter we always define k"^ > 0), which results in the appearance of the running 
coupling a(A;^). In an Abelian theory, due to the simple form of the Ward identities 
[i.e. the cancellation of the fermion wave function and vertex renormalization constants), 
the running coupling can be reconstructed using dispersion techniques (see [|18[)- The 
dispersion relation for the running coupling reads 



Jo jJi^ + k^ 



a(0) +k 



2 j--- dfj^ _p(jj^ ^^'"^^ 

/i^ /i^ + /c^ ' 
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with the photon wave function with the normahzation ^3(0) = 1. The spectral 
function p(/i^) (positive in QED) has support only on the positive real axis and is obtained 
by considering all discontinuities associated with the virtual photon. It is given by the 
discontinuity of ^3 at positive virtuality, 

Pif^')^^[Mf^'-^e)-2sif^' + ^e)] = -ioisc (2.2) 



The dispersive representation ( p.lD has the following non-perturbative implication: the 
coupling a{k'^) is regular in the full k^-complex plane with a cut on the negative real 
axis. Therefore this representation gives a coupling that is free from the spurious QED 
Landau singularity and deviates from the standard perturbative e.m. coupling a^'^{k'^) 
by a negligible amount O {P/A^) (A ~ lO^o-^OGeV). 



Non-Abelian case. Consider now the exchange of a gluon with negative virtuality —k"^, 
which should result in the appearance of the running coupling as{k'^)- In the QCD case one 
again seeks to identify the Feynman diagrams responsible for the running coupling with the 
one-gluon reducible graphs, which include gluon and quark self-energy insertions together 
with vertex corrections. Such a class of diagrams does not constitute a gauge-invariant 
set. Their gauge-dependent part is cancelled by corresponding contributions from the 
one-gluon irreducible graphs {e.g., two-gluon exchange). Therefore one can argue that 
these gauge-dependent contributions from the one-gluon reducible set of graphs remain 
finite at k"^ 0, while the gauge-invariant part behaves as contributing to the 

renormalization of the pole of the gluon propagator. Since we are going to concentrate on 
the leading behaviour at small k"^, we shall not analyse the contributions that are regular 
for fc2 _> 0. We emphasize again that in non-Abelian theories the very notion of one gluon 
exchange becomes elusive at the level of relative fc^ corrections, at which two gluons can 
mimic the one-gluon exchange, the two processes being comparable and related by gauge 
invariance. 

Inspired by the Abelian theory, we shall assume that in QCD the running coupling 
as(fc^) still satisfies a (formal) dispersion relation of the form (|2.1| ): 



00 , 2 1 
as{k') ^ aMZ{-k^) = -j Psif^') , Psip') = -— Disc {a,( V)} • (2.3) 

We note that the assumption of such a dispersion relation (with causal support) implies 
the absence of the perturbative Landau pole at k"^ = A^, the QCD scale, but it also 
involves such ill-defined quantities as the coupling ^^(O) and the spectral density ps at 
small scales /i^. However, as we shall see shortly, a<j(0) will not appear explicitly in 
physical observables, while the contribution of the spectral density in the small p"^ region 
will be suppressed by powers of /i^. 

Since in QCD the cancellation between the vertex and fermion propagator corrections 
due to Ward identities [Z^Zq = 1) does not hold, the quantity Z in ( p.3|) is not simply 
the gluon propagator correction. As is well known, both the non-Abelian part of the 
vertex renormalization correction Zj^^'' and the gluon propagator factor Zg contribute (in 
a gauge- dependent way) in forming the running coupling. Therefore in the case of QCD 
the quantity Z in ( p.3| ) is now given by the product 

Z = . Zg ■ 4°'^ . (2.4) 
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This expression shows that the function Ps(Ai^) is not hterally a spectral density (although 
we shall continue to call it by that name) since it is not necessarily positive. In QED 
the function p{fi^) is given by the discontinuity of the photon propagator, so that it is 
positive for any value of /i^ as required by unitarity. This implies that a(/c^) increases as k"^ 
increases. In QCD instead as(fc^) decreases with increasing k"^ and so Ps(/^^) is negative, at 
least in the perturbative domain. This is not in conflict with unitarity since, according to 
( ^^ , the function Ps(/^^) is given by the discontinuity of the gluon propagator correction 
(Zg) together with that of the non-Abelian vertex corrections (4"^^^). For example, in a 
physical gauge (axial, planar; see [1^) the first is determined by the total gluon splitting 
probability. 

Disc Z„oc dzl2Nc[- + 2 + z(l 

Jo [ \z 1-z / ^ 

and gives an infrared- divergent contribution, which is positive as required by unitarity. 
The second (vertex corrections) gives a (divergent) negative contribution, 

2 Disc 4°^^) oc -AN, ''''^^ 
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so that the total discontinuity is finite and gauge-independent but not positive definite: 

Disc [Zi^^^ ZgZi^^^) oc 1^ dzl 2N, (-2 + z{l -z)) + J2 {^^ + (1 - ^f) 

11 2 
= -—N, + -nf = -Po. 

2.2 Time-like gluon 

We consider now the case of a gluon with zero (fc^ = 0) or positive virtuality (fc^ > 0) 
corresponding to the contribution of new production channels. In inclusive quantities 
the integration over k"^ has to be performed up to some kinematic limit k'^^^ ~ Q^, the 
relevant hard scale. From this additional gluon correction one reconstructs the running 
coupling, provided one is able to factorize the following combination 

,,^o)sie) . ^ ^ ^iJ^] . (2.5) 

k'' TT [ —k'^ — te ) 

The first contribution is due the on-shell gluon while the second is due to the sum over all 
final states generated by a positive- virtuality gluon. In order to factorize this combination, 
we need to consider observables which are sufficiently inclusive that all new production 
channels contribute to the same value of the observable. 

Operationally our approach is similar to the so-called "naive non-Abelianization" pro- 
cedure based on resummation of fermion loop diagrams • The usual motivation for 
concentrating on fermion loop contributions is to avoid problems with gauge invariance. 
However we believe that the argument given earlier, that gauge-dependent terms are less 
singular at fc^ — > 0, justifies the application of our method to the leading power corrections 
in QCD. 

In the next subsection we describe how the running coupling emerges in inclusive 
quantities of this type. The dispersive representation ( |2.5| ) allows one to identify the 
proper scale for the argument of the running coupling, and to define a 'physical' scheme 
for the definition of the QCD scale A. 



2.3 Argument of the coupling 

For the sake of illustration we consider the non-singlet quark distribution F{Q'^, x) in DIS 
p!8| . Here p and q are the momenta of the incoming quark and of the hard photon probe 



F{Q\Qlx) = f^^kl l'^\e:F{Q\k\kl,x,z) L(0)5(A;2) + (2-6) 

Jx z JqI Jo V fc^ / 



respectively, so that = — and x = Q'^/2{pq) is the Bjorken variable. 

The evolution equation for F is obtained by considering the emission of an additional 
gluon of momentum k. For simplicity we discuss only the real emission contribution; 
virtual corrections can be included in a straightforward way. One has to take into account 
both the on-shell contribution at /c^ = 0, with coupling ^^(O), and the contribution from 
the continuum involving the discontinuity through multi-parton states, with strength 
given by the spectral density ps{k'^)- Introducing a minimum transverse momentum Qq 
for the emitted partons, the non-singlet quark distribution F{Q'^,Ql,x) is obtained by 
convoluting the combination (p.5|) with the matrix element squared T[Q'^ ^ k^,k'j_, x, z) for 
emission of an additional off-shell gluon of virtuality fc^, transverse momentum k^ and 
fraction of momentum {1 — z) = 2k ■ {q + xp) x/Q"^. We have 

'Ql 

where the integrations over k'j_ and the timelike gluon virtuality k'^ actually run up to the 
kinematical limit, W"^ = Q'^il — x) / x Q"^ . 

To see how the running coupling emerges in the anomalous dimension it suffices to 
consider the quasi-collinear region fc^ -C Q^. In the collinear limit the function can 
again be expressed in terms of the distribution F itself as follows 

The first term in the brackets, involving the standard q ^ q splitting function P{z) = 
{l + z'^)/{l — z), is related to the k^ = limit of the scattering matrix element. The second 
(subleading) contribution, non-singular for k"^ 0, originates from the /c^-corrections to 
the matrix element. 

In order to show how the integral over the virtual boson mass fc^ in ( p.6| ) reproduces 
an effective coupling running with fc^ as scale, it suffices to perform two steps, namely 
to extend the k"^ integration to infinity and then to apply the dispersion relation. The 
extension of the integration region is possible since the A;^ integral is rapidly convergent 
above k"^ ~ k'j_/z < k'j_/x which, in the collinear approximation, is much smaller than the 
kinematical boundary k'^ < ~ Q^. 

Given this simplification, one can use the dispersion relation ( |2.5| ) to obtain for the k'^ 
integrals in (|2] 



(2.J 



\kl + zk^ {kl + zk^yj y ' ^ ' ' A;2 

2 /l2 



asjkl/z) da,{kl/z) 
^ ' B dB ' 



and the evolution equation becomes 



Jx Z JQ^ K\ Z ZTl \ Z 



d\nk 



(2.9) 



Thus, by using the spectral function and the dispersion relation ( p.5D , we have recon- 
structed the running coupling associated with an extra gluon with positive virtuality. 
Moreover, as in |TP|,^,|22| , we have shown that the dispersion relation allows one to iden- 



tify the physical scale for the argument of the coupling. 

It is clear that in this way one has included some important contributions from higher 
orders proportional to the first beta function coefficient Pq. In particular in (|2.9|) one finds 
the following contributions: 



1. The last term in the curly bracket of (|2.9|) can be written as 

-'1 - -itir S"-'*^' ^ ■ ■ ■ • ("» = T^- 3"' 

2. Expanding the argument of the coupling in the main term, we obtain 

P{z) as{kl/z) = P{z) asikl) + P{z) In z (^^^ al{kl) + . . . 

The above terms of order are indeed present in the two-loop contribution to the anoma- 



lous dimension [23 



QCD scale. Two comments are in order concerning the "physical scale" of the coupling. 
1) Rewriting the expression in the curly brackets in (p.9|) as 



(2.10) 



one observes that beyond the first loop the coupling of the two pieces of the gluon radia- 
tion probability P{z) = 2z/ (1 — z) + (1 — z) effectively acquire different arguments. This 
is because these two contributions to the splitting function are physically different. The 
first term comes from the universal "soft" bremsstrahlung distribution {du/uj), which is 
independent of the nature of the incoming parton and corresponds to the gluon polar- 
ization in the scattering plane (longitudinal polarization), while the second one is due to 
"hard" gluons {u duj) , depends on initial parton spin and consists equally of longitudinal 
and transverse polarization. 

2) The universal nature of soft gluon bremsstrahlung may be used to define a "physi- 
cal" QCD coupling beyond the first loop. To do so one has to analyse the higher-order 
terms of the anomalous dimension and to absorb the singular contributions oc {1 — z)~^ 
systematically into a redefinition of the QCD scale A, which determines the intensity of 
radiation of relatively soft gluons. 

Such a physical scheme is related to the popular MS scheme as follows. In the MS 
scheme the intensity of soft bremsstrahlung is given by an infinite series in a formally 
defined parameter a^jg. In particular, both the first and the second loop MS splitting 
functions P{z) and P^{z) contain the "soft singularity" {l — z)~^. Redefining the coupling 
in such a way that the new two- loop splitting function P'^'^\z) does not contain the soft 
singularity (the so-called MC-scheme of [|^) one arrives, in second order, at the relation 



aj^P{x) + al^P^{x) = asP{x) + alP^'\x) + 0{o?^ ■ (2.11a) 
27r 



a, = 1 + ^ /C , A = A^exp(/C//3o) , (2.11b) 
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with 



which is equal to 4.0 ± 0.5 for = 4 =p 1. 

In this scheme the second loop contribution P^'^\z) is at least one power of (1— 2;) down 
with respect to the singular part of the main bremsstrahlung spectrum, P{z). Moreover, as 
advocated also in Ref. P^'^\z) no longer contains any n/-dependence, this ill-defined 
quantity being completely absorbed, to this order, into the momentum dependence of the 
coupling (for more details see p6|). 

Because of the classical character of soft bremsstrahlung, such a definition of the 
coupling proves to be universal with respect to the nature of the source of radiation. 
In particular the two-loop gluon splitting function also becomes infrared regular in the 
physical scheme described above, i.e., Pg'^gg{z) / Pg^gg{z) with z^l. For this reason 
we shall use this physical scheme for the QCD coupling. 



2.4 Effective coupling 

The spectral density ps may be related to the running coupling by the following formal 
transformation^ 

cysik') =-r -4^ Psif^') = - r 7^ Psivk^) 



/i^ + A;^ Jo 1 + V 

„ TT^-P '""^ ) (2-12) 



TT 



sin (tt^) 

In terms of the differential operator V, 



d_ 
(iln/i^ 



^^TfT-^' (2-13) 



one may thus write the inverse relation as 



p,(/i2) = -sin(7rP)a,(^2) _ ^2.14) 

77 



Introducing the effective coupling aefr(yU^) defined by 

d 



din p? 
we have 



Ps(p) = Ti ^ «efr(/i ) , (2.15) 



a.(/.^) = A:^/f P^M^^) (2.16) 



•^Let us recall that we are considering here the coupling constant Q!s(fc^) which is free from spurious 
singularities in the Euclidean domain, fc^ > 0. 



Q 



with inverse 

aesifi^) = ^ asifi^) • (2.17) 

It follows from Eq. ( p. 17 ) that in the perturbative domain <C 1, the standard and 
effective couplings are approximately the same: 

a^sift^) = _ + . . . = + O fa^\ . (2.18) 

b dm fi'' ^ ^ 

In what follows we shall look upon aes itself as an effective measure of QCD interaction 
strength which extends the physical perturbative coupling down to the non-perturbative 
domain. 



3 Dispersive method 

Consider a dimensionless inclusive quantity F{Q'^, {x}) for a hard process involving only 
quarks at the Born level. Here is the hard scale and {x} stands for any further 
relevant dimensionless parameters, e.g., Bjorken x or its conjugate moment variable N 
(DIS structure functions), particle energy fraction (inclusive e^e~ annihilation spectra), 
the ratio M'^ / s (the Drell-Yan process), or a jet shape variable. The one-loop prediction 
for Fi^Q"^, {x}) , obtained from the squared amplitudes involving one additional gluon, 
can be improved by using a dispersive method based on the representation of the running 
coupling in for the space-like case and ( |2.5| ) for the time-like case. In this way one 
takes into account relevant higher-order perturbative effects and is also able to study 
power contributions arising from the non-perturbative behaviour of the coupling at low 
scales. 



3.1 Basic equation 

To describe the method we concentrate first on collinear safe quantities. The generaliza- 
tion to include collinear singular observables will be presented shortly. We have to consider 
the two cases in which the additional gluon contributes to real production channels and 
to a virtual correction. 

Consider first the case in which the additional gluon contributes to new production 
channels. As discussed in Sect. 2.2, one has to take into account both the on-shell contri- 
bution at k"^ = 0, with coupling as{0), and the contribution from the continuum involving 
the discontinuity through multi-parton states (from gluon branching), with strength given 
by the spectral density ps{k'^). For all these contributions at fixed k'^ one would like to 
factorize the combination ( |2.5| ). This is possible only if the observable is fully inclusive 
with respect to gluon branching, as in the case of total cross sections, DIS structure func- 
tions, etc. For less inclusive quantities, such as jet shape observables, gluon branching 
may give a different contribution The non-factorizable contribution is higher-order 
in ctefT- However, since a^^ enters at a low scale, such terms could still be comparable 
to the factorizable first-order contribution oc Ues- Therefore our method can only give 
quantitative predictions for less inclusive quantities if at low scales is sufficiently small. 

We denote by T^'^^^{Q'^, {x}, /c^) the squared amplitude for the emission of a gluon with 
timelike virtuality fc^ > 0, apart from the coupling. JF^'^'^^(Q^, {x}, fc^) is defined to vanish 



o 



outside the real phase space < fc^ < fc^axlQ^' {^D- Then the real emission contribution 
to the inclusive quantity F{Q^, {^}) is given by the sum of the on-shell and continuum 
contributions (see Eq. (^.5D): 

r^^\Q\ {x}) = «,(o)^'^<=^'(o) + / ^psik'') ■ r'^\k^) . (3.1) 

JO K 

Consider now the case in which the additional gluon contributes to virtual corrections. 
This contribution is present when the observable under consideration does not vanish 
in the Born approximation. In this contribution the virtuality of the gluon is spacelike, 
— fc^ < 0, and the coupling is given by as(fc^). The contribution can then be written as 

dk 



with M. the integrand for the Feynman diagram containing the (massless) virtual gluon, 
apart from the coupling. We now use the dispersive representation ( p.3| ) for the running 
coupling to write 



virt /n\ / l2\ , "A'' / 2\ / "'^ \ /I t 1 2\ 



= a.(0)^-"*(0) + /°° -^p.in^) ■ r"'{p^) 

Jo /i^ 

where 



(3.2) 



^-^HQ^{x},/x2)= / ^^A<(g^{x},-P) (3.3) 
JO k-^ + /x^ 

is the one-loop virtual correction that would be produced by a gluon with a finite mass 

The result thus has the same structure both for the real (|3.1|) and for the virtual (p.2|) 
contributions. Combining them to form 

which we shall refer to as the characteristic function, the improved one-loop formula reads 
(omitting the external variables and {x}) 



F = as{o)m + r ^Psif^') ■ Hf^') = r %Ps{p') ■ [H^^') - m ] , (3.4) 

JO p^ Jo p'^ ^ J 



where we have made use of the formal relation (|2.3| ) to eliminate q;s(0). By introducing 
the effective coupling aeff(/i^) ( |2.17] ) using Eq. ( |2.15| ) and integrating by parts, we can 
write 

F{Q\ {x}) = r ^ aeff(/i') ■ HQ\ {^}; , ■ (3.5) 

JO p'^ a m /i^ 



The relation ( |3.4] ) or ( |3.5| ) is our basic equation for studying both the perturbative part 
and non-perturbative contribution to F. For a coUinear safe observable ^{p"^) vanishes for 
/i^ — > and /i^ — > cxo, and the integral is dominated by the region p^ ^ Q'^. One observes 
that the integrand contains two well separated scales: Q^, the scale of the characteristic 
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function J-'{fi'^), and A^, the QCD parameter, which is the scale of the effective couphng 
Q^efr(/^^) or Psifi'^)- The perturbative contribution is related to the scale Q^, while the 
power-behaved contribution is related to A^. 

Let us stress that although in our derivation we have introduced the quantity as{0), 
which may make little sense in the context of QCD, the final results (|3.4|) or (|3.5|) involve 



only the spectral density Ps(/i^) or the effective coupling aeff(/i^) convoluted with a func- 
tion that vanishes at /x^ = 0. Therefore the basic equation may make sense even if as{0) 
does not exist. 

In what follows we shall look upon a^s defined in (|2.171 ) as an effective measure of 
QCD interaction strength which extends the physical perturbative coupling down to the 
non-perturbative domain. While the form of aes{k'^) at large k"^ is well known, we treat 
OefT at small scales as a phenomenological function whose behaviour is to be determined 
from experiment. 

The characteristic function is process-dependent and is obtained in practice by 
computing the relevant one-loop graphs with a non-zero gluon mass /i. Note that we do 
not intend to imply that the gluon has in any sense a real effective mass, but only that 
the dispersive representation can be expressed in this way. 

In the above derivation it was implied that the characteristic function is well-defined at 
/i^ = 0. This is true for coUinear safe quantities such as the total e^e~ annihilation cross 
section, DIS sum rules, the Drell-Yan K factor, jet shape variables, etc. On the other 
hand for coUinear singular quantities, such as DIS structure functions, the Drell-Yan cross 
section, and inclusive e^e~ parton fragmentation functions, we have 
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^(g^ {x}; = P{{x}) ■ In % + r^>^{Q\ {x}; /i^) , (3.6) 



where J^^^^{Q'^, {x}; 0) is finite. We cannot use the representations (|3.4|) for the collinear 
singular distribution itself, but rather for the scaling violation rate Q^dF / dQ^ . In this 
case the function is involved, which has a finite /i^ — limit. Then Eq. ( p.5|) becomes 

where, from ( |3.6| ) we have 

TiQMxhf^') = (^6^^ -^^^^(e, {x}) , (3.8) 

and we have used the fact that JF is a function of e = fi'^/Q'^. The function JF vanishes as 
e ^ and e —>■ oo, and therefore the integral in (|3.7|) is now well defined. For simplicity, 
we continue the discussion in terms of JF and JF, with the understanding that one should 
substitute JF and JF when studying the scaling violation in a collinear singular quantity. 

We now discuss the general behaviour of the characteristic function. Since JF depends 
only on dimensionless ratios, we write 



^(Q^M;/i2)=^(e,{x}), ^^-e^^(e,{x}), e^^. (3.9) 



We describe the form of the characteristic function JF(e, {x}) for a collinear safe quantity, 
and in particular its behaviour for large and small e. For a collinear singular quantity, 
the same behaviour is obtained for its regular part jF'''^s(e, {x}). 
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Large e. At large e the characteristic function either vanishes identically because of the 
phase space boundary (in the absence of virtual corrections to this order, as is the case 
for jet shape variables), or decreases as a negative power of e as a consequence of the 
renormalizability of the theory. It follows that the logarithmic derivative T also vanishes 
at infinity, so that the integral in Eq. (|3.5D is well-defined. 



Small e. The behaviour at small e is crucial for the analysis of non-perturbative power 
contributions. By definition of a collinear safe quantity, the behaviour of T near e = is 
of the form 

j={t) = [f{\ne)+eg{\ne) + 0{e')] , (3.10) 
with J9 > 0, where / and g are polynomials of degree not higher than 2: 

/(lne) = /2ln2e + /ilne + /o, (3.11) 

and similarly for g{lne). 

As an example, we show in Figs. and ^ the behaviour of 7?. = 2-71 J-" /Cp and its 
derivative TZ for the e"*"e~ total annihilation cross section. In this case we have TZ ~ 
3e^ + 2e^ In e — 3e^ at small e, as shown by the dashed curve in Fig. 0, so p = 2, /(In e) = 3 
and g(\n e) = 2 In e — 3. 

We show finally how the standard one-loop result is recovered from Eq. ( p.5|) . Since 
JF vanishes both for e and for e ^ oo, the main contribution to the integral comes 
from e ~ 1, that is yU^ ~ Q^. Therefore one can extract the leading contribution by taking 
the value of the effective coupling at some characteristic scale 

^ c({x}) ■ 

outside the integration, to obtain 

F(Q^ {x}) = aeff(Q') ■ ^(Q', {x}; 0) + T ^ [ aefr(/i') - «efr(g') 1 ■ HQ^ i^}; /^') • 

Jo fl^ ^ J 

(3.12) 

The remaining integral contributes both to higher-order perturbative terms and to non- 
perturbative power contributions. Roughly speaking, the former come from the region 
/i^ ~ (Q^ ~ and the latter from /i^ ~ A^ ^ Q^. The small- /i^ region is outside 
perturbative control but, to the extent that aefr may be defined universally and obtained 
from experimental data, its contribution is determined by the behaviour of JF at small 
values of e. In the remainder of this Section, we consider various approaches to the 
evaluation of the integral in ( p.l2| ). 



3.2 Renormalons 

For large the main contribution to the integral in (p.l2| ) still comes from the region 
/i^ ~ Q^. One might take this as sufficient motivation for evaluating this contribution by 
expanding aesifJ^"^) around /i^ = (for simplicity we consider here = Q^)- In this 
region the effective coupling aes{^'^) can be approximated by its one-loop (or two-loop) 
perturbative expression a^^^fi"^) 
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Making the replacement 



aefr(/i^) a^i^^') = + £ f^^^ V ' (3-14) 



k=l 



one finds 




(3.15) 



where the simphfied notation T{Q \ {x}, /i ) = Tie) has been used. The deficiency of the 
perturbative expansion becomes apparent from the fact that the resulting coefficients 



exhibit factorial growth . This is associated with both the ultraviolet and the infrared 
integration regions in e. 

The ultraviolet contribution to is estimated by integrating ( p.l5| ) over e > 1. If 
JF(e) vanishes like at large e, one finds for large k 

Cuv^ A-rfe / ly ^ (3.16) 
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This corresponds to an ultraviolet renormalon. Such an alternating series can be evaluated 
by Borel summation. This is because in the ultraviolet integration region of ( |3.12| ) the 



replacement aefr(/i^) — *■ oF^{\x^) is a reliable approximation and the contribution of this 
region can in fact be evaluated explicitly without any expansion. 

The infrared contribution to is estimated by integrating over the region e < 1. 
Using the small e behaviour in (p.lO ) one finds 



C^f = /^(lni)' 6^/(1116) ~ k\ (3.17) 

One again finds a factorial behaviour (an infrared renormalon). In this case however the 
coefficients are non- alternating and therefore the series is not Borel-summable. Attempts 
to ascribe meaning to such a series by brute force tend to give rise to unphysical complex 
contributions at the level of Q"^^ terms. This is generally interpreted as an intrinsic 
uncertainty in the summation of the perturbative seriesQ. It is important to recognize 
that formal mathematical manipulations alone cannot resolve this problem, which is of 
a physical nature. One requires genuinely new physical input, namely a power-behaved 
'confinement' contribution, to obtain a sensible answer. In this paper we advocate the 
hypothesis that such an input, for sufficiently inclusive Minkowskian observables, may be 
embodied in the form of aeff(/i^) at small /x^. 

3.3 Soft confinement and power corrections 

The standard operator product expansion (OPE) approach by Shifman, Vainshtein and 
Zakharov [|l^ quantifies confinement effects in terms of additive contributions to Euclidean 

"^In fact, infrared renormalons are a purely perturbative phenomenon and have no direct relation to 
the presence of the Landau singularity in the running coupling; for a detailed discussion see 129]. 
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quantities and supplies the fuel for the impressive machinery of the ITEP sum rules (for 
a review see [^] ). 

The OPE approach is based on the hypothesis that the entire effect of confinement 
in vacuum correlators of currents may be embodied into "condensates" (vacuum expec- 
tation values of gauge- and Lorentz-invariant colourless operators built from gluon and 
quark fields). The basic ITEP idea was to separate the long- and short-distance contri- 
butions and to treat them on different bases. In the Euclidean region such a separation is 
straightforward: for A;^ = /cq -|- > ~ 1 GeV^ one employs perturbation theory — in 
particular, one uses the purely perturbative expression for the coupling in Feynman 
diagrams for a given correlator. The complementary region k"^ < A^, on the other hand, 
is treated phenomenologically. Upon integration, the latter region gives a power-behaved 
contribution of the order of (A^/Q^)^- 

As far as gluon propagation is concerned, within the logic of the present paper one 
may write equivalently 

a.iP) = ank^) + 6as{k^) , (3.18) 

with Sag a modification in the effective interaction strength at small momentum scales 
responsible for non-perturbative effects. In agreement with the ITEP point of view, while 
this contribution generates power corrections to the hard distributions, it should not 
modify the behaviour of the running coupling at large momentum scales. 

Let us stress that the separation ( p.l8| ) does not mean that the perturbative contribu- 
tion as such is free from power-behaved terms. This very question is almost meaningless: 
one would have to keep under control an infinite series of logarithmic high-order terms 
prior to addressing the problem of perturbatively-generated power corrections. As we 
have mentioned above, due to the "infrared renomalon problem" the perturbative series 
is apparently not summable to such a level of accuracy. At the same time, the separation 
prescription (p.18 ) cures the problem operationally: as long as perturbative integrals are 



cut off in the infrared, resummation of the infinite series triggered by the running of a^"^ 
is harmless. To obtain a physically sensible (and reasonably accurate) answer it suffices 
to calculate a few terms of the perturbative expansion and then to add power-behaved 
terms as a contribution of essentially different (non-perturbative, confinement) origin. 

There is a strong point of the OPE ideology, which one may refer to as the "so/i 
confinement" scenario. The OPE prescription implies that the propagation of quarks and 
gluons with large (Euclidean) momenta remains unaffected by non-perturbative physics, 
even at the level of power-suppressed terms. Formally, one might argue that the propaga- 
tors of coloured fields have no gauge-invariant meaning. Nonetheless a power-suppressed 
variation, say a 1 + A^/fc^ correction to the gluon propagator in the ultraviolet region 
k^ ^ X^, would inevitably introduce an additional ^^non-perturbative" contribution. This 
would have no relation to the region of small momentum flow in the corresponding Feyn- 
man diagrams for the current correlators, which region is responsible for the formation 
of condensates. As a result, the original separation idea would not be valid. Thus the 
term 5as(^^) in (|3.18|) should decay faster than some (sufficiently large negative) power 
fQj^ viability of the notion of a condensate of dimension 2p < 2pmax- 
Physically the OPE prescription corresponds to a picture of smooth non-perturbative 
large-scale vacuum fields with a typical size ~ A^^. If such fields were non-singular at 
short distances, the gluon propagator (and thus the running coupling) would be subject to 
exponentially small corrections only. An instanton-anti-instanton g 3jS clS 8b representative 
model for non-trivial vacuum fields sets an upper bound Pmax < /?o ~ 9, above which 
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small-size instantons start to disturb the propagation of quarks and gluons. 

Following the discussion in the previous sections, we are now in a position to implement 
a similar logic for quantifying confinement contributions to Minkowskian observables. 

The effective coupling aesifi"^) appearing in 
turbatively below /i^ 



and (|3.5| ) cannot be defined per- 
1 GeV^, where the very language of quarks and gluons is 
scarcely applicable. In this region we expect an "effective coupling modification" 6aes 



which generates the non-perturbative interaction strength 6as in (|3.18| ) via the dispersion 
relation ( p.l6| ), i.e. 

Sasi^) = e r , 5aefr(/i') • (3.19) 



h (/i2 + P)2 

At first sight, one might feel free to choose an arbitrary form for 5aefr at small scales 
to model confinement effects. However, on inspecting ( p.l9| ) one observes that, generally 
speaking, a finite modification of the effective coupling at low scales will affect the ultra- 
violet behaviour of the coupling in the Euclidean region by an amount proportional to 

As discussed above, such a modification would ruin the basis of the OPE approach. 
One has therefore to require that at least the first pmax integer moments of this coupling 



vanish. A similar constraint has been discussed in 30 



Consider, for example, the non-perturbative gluon condensate which contributes to 
the Adler D-function (see, e.g., [0). To first order in it is given by the integral 



TT / 



NP 



3CV 



Jo 



(3.20) 



where is the ultraviolet cutoff. Substituting the representation ( p.l9|) for 6as in terms 
of the function Saes and performing the integration over k"^, one obtains 



27r2 
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NP 



3Cp 
Jo 



Convergence of the integral in ( |3.20|) implies 



{f^'Y Sa^^if^') = 0; p=l,2, 



/o /i 

that is, the vanishing of the first two moments of 6aes{f^'^)- The result now reads 



(3.21) 



TT 



9 \ TT /NP 



2tc Jo li^ 



dfi 



(3.22) 



Notice this (log-)moment integral is independent of the scale of the logarithm. 

Inspecting higher moments (operators of higher dimension), one finds it necessary to 
impose on Saesi^fj,"^) the set of restrictions ( |3.21| ) with p = 1,2, . . .pmax, to respect the cri- 
terion of "soft confinement" as seen through Euclidean eyes. Applied to a Minkowskian 
observable F, this means that only those terms in the small-yU^ behaviour of the character- 
istic function J^(/i^) that are non- analytic in /i^ will contribute to F^^. Among them are 
the terms with non- integer p and/or those with log /i^-enhanced asymptotic behaviour. 

Eq. (|3.22| ) gives an example of a "measurement" of one of the log- moments of 6aes. 
In the rest of the present paper we study how this and other non-zero moments of 6aeff 
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enter into different Minkowskian observables. According to Eq. ( |3.5| ), the corresponding 
contribution to a generic collinear safe observable F will be of the form 

FNP(Q2, {x}) = Sa^sifi') HQ', {^y, /i') • (3.23) 

The leading non-zero contribution to the integral will come from the first non-analytic 
term in the small-e expansion of J-'{Q^,{x}; fi"^) = jF(e, {x}), Eq. ( |3.10| ). If p is not an 
integer, the first term is non-analytic and the leading contribution is proportional to 
Q-2p It^q q -^j^ere fq In'^ e is the largest non- vanishing term in Eq. ( ^.llj ). If p is an integer 
and /i or /2 is non-vanishing, the leading contribution is proportional to Q^^^ln'^^^ Q, 
because the first singular contribution in fi"^ will be of the form /i^^ln /i. If p is an integer 
and both fi and /2 are zero, the contribution will be proportional to In"^"^ Q, 

where g > 1 is now specified by the largest non- vanishing log-enhanced term of g{t). 

We shall call power-behaved contributions of this type dispersively- generated power 
corrections. An important question to be addressed is: how much of the full power 
correction to a given observable is due to these dispersively-generated terms? 

For quantities like DIS structure functions and the Drell-Yan i^-factor, dispersively- 
generated l/Q^ power contributions arise which are given by the first (log- and log^- 
enhanced) moment of Sacsif^"^)- The origin of these contributions may be traced back 
to a universal l/Q^-suppressed non-perturbative correction to the relevant hard cross 
sections (coefficient functions). At the same time, from general OPE considerations one 
expects l/Q"^ corrections proportional to the hadronic matrix elements of the relevant 
twist-4 operators, which depend on the target hadron. Since the dispersively-generated 
terms have a definite, calculable dependence on the hard process kinematics, it may be 
possible to identify kinematic regions in which they are dominant. 

In the case of e~^e~ -annihilation, which is a hard process free from initial-state hadrons, 
the situation might be simpler. In particular, for the practically important case of linear 
(l/Q) corrections to event shapes, dispersively-generated terms could describe the entire 
contributions to different quantities in a universal way. Such a hypothesis seems plausible 
not only because there is no competition from the conventional OPE contributions, but 
also because it follows naturally from the picture of soft (local in phase space) hadroniza- 
tion, which one invokes to explain the observed similarity (duality) between the calculable 
distributions of partons and the measured distributions of hadrons from QCD jets ||3T[|. 



4 Applications 

In this section we apply the dispersive method of the previous section to various hard 
processes dominated by quarks. Recall that the object of central importance is the char- 
acteristic function J^{e) for the emission of a gluon with mass-squared = eQ^ at the 
hard scale Q^. For power corrections, the relevant contribution is given by the leading 
non-analytic term in the small-e behaviour of the logarithmic derivative JF, which is of 
the general form ( p.lOp . The coefficient of the resulting contribution is then determined 
by the function /(Ine) (or by the next-to-leading function 5'(lne) if the leading term is 
analytic). Since / or (7 is a polynomial in ln{fi'^/Q'^) of degree not more than two, we 
introduce for future use the moment integrals 

= ^ r ^ /i^^ 5aefr(/x^) (4.1) 
Ztt Jo u'^ 
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and their derivatives 



A" - A.. 



(4.2) 



which have respectively an extra factor of In and In^ /j^ in the integrand of Eq. ( [4. 1|) . 
Then all dispersively-generated power corrections proportional to Q~^^ can be represented 
as linear combinations of A2p, A^p and Ag^. Since the integrand must be non-analytic in 
/x^, the unprimed moments can only contribute when p < Pmax is not an integer, while 
the primed (log) moments can contribute for any p. The quantity on the right-hand side 
of Eq.( |3.22| ), for example, is denoted by A'^, corresponding to a l/Q"^ correction to the 
Adler D-function. 

In general the characteristic function J^{e) has both real and virtual contributions. 
Since the latter are universal we discuss them first. Their form depends on whether the 
momentum transfer in the hard process is space-like or time-like. 

Consider the one-gluon virtual correction to a hard process {e.g. DIS) with a space- 
like momentum transfer —q^ = > 0. The total correction to the renormalized hard 
interaction vertex is of the form CFasVs/^ir where 



v.. 



de g2 



-2 dz(l - zf (;.2 + ^2)(p + ^(J2) 



_2 r,z^i^i.i 

Jo z — e e 



2(1 -e)^ 
2(1 -e)^ 



Li2(l - 1/e) 



TT 

~6 



7 



(3-2e)lne + 2e 



1 vr^ 

Li2(e)+lneln(l-e)--ln\-y 



(4.3a) 



^ - (3-2e)lne + 2e; 



where 



U2iu) = - r^ln(l-t). 

JO t 



For a time-like process (g^ = > 0, as in e~^e annihilation and the Drell-Yan process), 
the virtual correction is Ci?asVt/27r where 



Vtit) = Re Vs(-e) = -2 / dz^ ^ In - 

JO z + e e 

2(1 + e)^ Li2(-e) + In e ln(l + e) - ^ In^ e + ^ 

2 



(3 + 2e) lne-2e, 



(4.3b) 



Next we combine these with the real contributions for various hard processes and analyse 
the resulting power corrections. We consider first collinear safe and then coUinear singular 
processes. 

4.1 Total e+e" annihilation cross section 

The most straightforward application of the method is to calculate the power correction 
to the e^e~ annihilation cross section R{Q'^), given to first order by (we suppress the 
standard parton model normalization factor = NcJ^j^}) 



+ ... 



(4.4) 
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In our approach, this result originates from perturbative evaluation of the expression 



R 



1 + ^ / ^«efr(/i')■7^(Q^/i2^ 



n 



d 



din /i' 



(4.5) 



where TZ is the characteristic function obtained from diagrams involving one gluon with 
mass yU. This function has both a real and a virtual contribution (see Fig. ||). The latter is 
simply given by the function Vt(e) in ( 4.3b ). To obtain the real contribution we consider 
the emission of a quark, an antiquark and a gluon of momenta p, p and k respectively 
(fc^ = /i^, = = Q). The matrix element assumes a simple form in terms of the scaled 
quark and antiquark energies, 

X = 2pQ/Q^ , X = 2pQ/Q^ , 

which satisfy the phase-space constraints 

(1 — x)(l — x)>e, X + X < 1 — e . 

The matrix element squared is then of the form CpasAiee/^TT, where 

e e 



MeJx,x,e) 



(x + eY + (x + e)2 



• ■ (l-x)(l-x) {1-xy (1-x) 
The real contribution to the characteristic function TZ is thus 

dx dxA4ee{x, X, e) 



(4.6) 



= -2(1 + 6 



ph.sp 
2 



2Li2(-e) + 2 In e ln(l + e) - - In^ e + ^ 

2 



+ 5 + (3 + 4e + 3e^)lne-5e2, 

(4.7) 



with e < 1. At the edge of phase space (e 1) the distribution vanishes rapidly, as 

7^«(6) = l(l-e)^{l + 0(l-e)}. 

Finally the complete characteristic function is 

for e>l, 7^(e) = Vt(e), 
for e< 1, 7^(e) =7^('')(e) + V^(e) 

= -2{l + ef[ Li2(-e) + In e ln(l + e) ] + - + (2 + 3e)e In e - 2e - Se^ . 



(4.8) 



(4.9) 



The function TZ{e) is plotted in Fig. |T| together with its logarithmic derivative TZ in 
Fig. 0. The behaviour for small and large e is as follows: 
For e < 1 
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+ O(eMne) 



7^=:^(lne + ^) + O (^-^ 
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(4.10a) 



(4.10b) 



These two limiting forms are in fact related by the following symmetry of the e+e 
characteristic function. 
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Figure 2: Derivative of characteristic function for e'^e~ total cross section, TZ. Dashed 
and dot-dashed curves show the hmiting behaviour at small and large e, respectively. 



1 n 



Symmetry. The function TZ given in ( [4.9| ) satisfies the following inversion symmetry 



TZie) 



7^(l 



(4.11) 
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The behaviour for e — > follows from this symmetry and the fact that for e 
characteristic function 7^(e) = Vt(e) vanishes in the way specified in ([4.10b| ). 

The symmetry ( |4.11| ) also explains the fact that in the sum of real and virtual contri- 
butions one finds that not only the terms singular as fi"^ —>■ 0, 



InV 



In/i^ 



cancel, as required by the Bloch-Nordsieck theorem, but that the cancellation extends 
||7|JT7[| also to the following finite terms 



fx^ln fi , /i^ In /i . 



/i^ In /i , /i^ In /i . 



Perturbative evaluation of R{Q^) . Now we are in a position to evaluate R according 
to ([4.5|). Since TZ vanishes both for small and large e. 



7^ = Se^ + 2e^ In e + . . . , e 



7^ 



3e 



In e + . . . e ^ oo , 



(4.12) 



the main contribution to (^4.9|) comes from e ~ 1, that is /i^ ~ (see Fig. |^). Therefore 
the leading perturbative contribution to R^^{Q^) is obtained by taking the value of the 
coupling constant acsiQ"^) outside the integration and one finds 



1 



efr(Q')7^(0) + 0{a',{Q' 



(4.13) 



where TZ{0) =3/2 and ~ should be chosen in the vicinity of the peak, so that 
c^esiQ^) ~ CisiQ"^) and one obtains the result ( f4.4|) . 



Leading power correction. As indicated in Eq. ([4.12|) , the first non-analytic term 
in the expansion of TZ at small e is of order e^lne. Thus the leading power-behaved 
contribution is given in terms of the moment integral (^4.1|) for the effective coupling as 



^6 

2— . 



(4.14) 



The absence of a l/Q^ contribution follows from the lack of any suitable dimension-two 
operators, in the massless quark limit. In principle one might have expected that al/Q^ 
contribution could be present, due to the gluon condensate (a^G^). As we have already 
remarked, such a contribution is indeed expected in the Adler D-function, but to first 
order in 6aes it does not appear in R itself, which is related to the discontinuity of the 
D-function In order for a.l/Q'^ term to appear in R one would need a In Q^/Q^ term 
in D. At the same time, a leading power correction of the form ( [4.14| ) is consistent with 
the OPE, since the Z)-function does acquire a log-enhanced contribution from operators 
of dimension six (see IIT^JT3 



on 



Hadronic width of the r lepton. The r decay width [p|-[T^ is closely related to the 
e+e" annihilation cross section. It is normally expressed in terms of the quantity Rr, 
which for massless quarks is 



Rr 



-I ds 
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Thus we can write 



where 



TZr {ml, /i^^ 



1 + 



2n Jo /i^ 

? ds 



mt 



mt 



mt 



Defining y = fi'^/m'^ and e = yU^/s, we can write this in the form 



T^Ay) = 2y 



de f^_y^'^ 



1 + 2 
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e J 



] 7^(e) . 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



The form of TZriy) is shown by the solid curve in Fig. ^ At large y the behaviour is 

'Iny' 



y 



shown by the dot-dashed curve. At small y 



n^y) = 8 (4 - 3C(3)) y - W + 6 In^ y - 12 Iny + - h/^ + 



32 



(4.19) 



(4.20) 



(dashed). 

The fact that the non-analyticity of the expansion ([4.20|) starts at order y^ is again 
in accord with the absence of a dimension-four contribution in the operator product 
expansion ||l2| , pr5| . In this case the In^ y factor indicates a logarithmic enhancement of the 
power-behaved correction to Rr-, which takes the form 



R. 



NP 



-6 (2A; In ml + 2^^ - A!^)/ml 



(4.21) 



4.2 DIS structure function F2 

We next consider the deep inelastic scattering process in which a quark, with momentum 
p, is probed by a hard spacelike photon of momentum q (g^ = — Q^, x = Q'^/2pq). Recall 
from Sect. 3 that in the case of coUinear singular quantities, such as DIS structure func- 
tions (or their moments) and e~^e~ fragmentation functions, the characteristic function 
J-'{Q'^ , {x}; ^'^) has a finite Q^-independent coUinear limit, J-'{Q^, {x}]0) = P{x). For 
small e (inside kinematical limits), 

^(e) =P(x)+ j-,eg(e), (4.22) 

where J-'reg{^) vanishes as a power of e at e — > 0, as does J-' for the coUinear safe case, see 
( p.lOp . As a consequence, the function JF vanishes as a power for both e — > and e — > 00. 




This quantity determines the scahng violation according to Eq. ( |3.7| ). More precisely, 
defining the moments of the non-singlet part of the structure function F2 as 

F2(Q2, N) = f' x''-'F2{Q\ x) dx , (4.23) 
Jo 

the scaling violation is described by 

T^{Q')^Q'j^\nF,{Q\N) = ^ T ^ aeff(/i') J^(e = /iVg^ iV) . (4.24) 
oQ"^ Ztt Jo jj,'^ 

Evaluation of characteristic function. The characteristic function T{e) is obtained 
by considering the amplitude for an incoming and outgoing quark and for an off-shell 
gluon of momentum fc, either virtual or real. For the combination of polarizations leading 
to the F2 structure function, the amplitude-squared for off-shell gluon emission assumes 
the following simple form (apart from the coupling as and the quark colour factor Cf/'^t^) 

y(l — x) — xe 2a;(l — y)(l — e) y(l — x) — xe ^ y'^(l — y) , . 

M2 = + ( + \2 + 4 ' (4-25) 

[y — xe)"^ [y — xe)[l — X) [1 — x)"^ [y — xey 



with y = {pk)/{pq). The last term in ( |4.25| ) is the one contributing to the longitudinal 



structure function F^ ■ When e = the first and second terms are collinear singular for 
y^O, with coefficient equal to the quark splitting function P{x) = (1 + a;^)/(l — x). 

The real part of the characteristic function contributing to F2 is obtained by integrating 
y over the phase space region ex/{l — x) < y < 1. Including the virtual contribution, 
which is now Vs(e) given by Eq. ([4.3|), one finds that 



J^(x; e) = J^(")(x; e) 6(1 - x - ex) + V,(e) 5(1 - x) (4.26) 



00 



with 



2(1 



- (1 + x) + 2(2 + a; + Qx')e - 2(1 + x + 9x^)e^ 



1 — X 

3 + 14e - 15e2 



In 



;i-ex)(l-x) 



+ 



+ 



+ 



X 



2(1 -x) ■ (l-x)2 2(1 -x)3 ■ 1-ex 
+ 1 + 3x + 6(1 - x)(l + 3x)e - (8 + 9x + 18x^)6^ . 

(4.27) 

The coefficient of — Ine is the quark sphtting function P(x), which is singular for x — > 
1. This singularity is regularized by including the virtual contribution. The F2 quark 
structure function to this order is thus (we suppress the standard parton model factor 



5(1 -x) + — 



Cf r°° dfx^ 



(4.28) 



2tt Jo /i^ 

The integral is convergent for fi"^ 00. In this limit the characteristic function J-'{e) is 
given by the virtual contribution alone and we have 



2 /, 11\ /lne\ 



3e 



6 



(4.29) 



In Figs. ^ and |^ we show the ffist few moments of JF and JF, for = 2, 3, 4. Note that 
the N = 1 moment of JF vanishes owing to the following identity, which holds for any e, 



Jo 



(4.30) 



This means that the Adler sum rule is satisfied identically, that is, it receives neither 
perturbative nor power corrections within our approach (see 



Logarithmic scaling violation. The moments of the limiting function P(x) in (|4.22|) 
have the meaning of the anomalous dimensions 'jN-f\ The usual perturbative scaling 
violation result for the moments of the DIS structure function, namely. 



Cat (a J 



is then reproduced as follows. Since the integrand in 
one obtains the leading contribution by substituting aefr(/^^) ctsiQ^)'- 



(4.31) 



is peaked at //^ = ~ 



^(1) 

N 



Zn Jo yU 

c^siQ')^ \HQ', N; 0) - /(Q^ iV; 00) 1 = a.m^PN 

/TT L J ZTT 



(4.32) 



where P/v is the corresponding moment of P(x). Keeping the next term in the expansion 
for the running coupling. 



^It should be noted that in higher orders P{x) will start to depend on the integration scale through 
the effective coupling, P = a^ff (/i^)), while remaining Q^-independent. 
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OA 



the first correction to ( [4.32| ) can be derived as 



where 

= lim |^(Q^ iV; _ In ^ . I = J-,,3(Q^ iV; 0) (4.34) 

is the coefficient function, the value of which clearly depends on the choice of the expansion 
scale Q^. An obvious (ffist-order) identification is as follows, 

7Ar(as) asCpPN/'^T!' , Ciyias) <^=^ asCpCN/^n . (4.35) 

Notice that in the original expression ( [4.24|) there is no arbitrariness in the choice of the 
hard scale (scheme dependence), which only emerges when one tries to evaluate the area 
under the JF curves in Fig. |^ (weighted with aesifi^)) in terms of a series expansion around 
a given point Q^. 

Another point that is clearly seen from Fig. ^ concerns the actual hardness scale of the 
process. The region from which the scaling violation rate receives its main contribution 
shifts to lower momentum scales with increasing N; in fact the peak is around ~ Q'^/N, 
or, equivalently, (5^(1 — a;) ~ W^. To avoid confusion, let us stress that identification of 
W'^ with the proper physical hardness scale is true for non-singlet structure functions 
only. Therefore it should not be applied to the small-x region, which is dominated by 
the singlet contribution, the scale of which does exceed but by far less than the 
kinematically allowed limit W"^ ^ Q'^/x ^ Q^. 

To obtain the power corrections we have to consider the small-e limit of the charac- 
teristic function. 



Small-e behaviour. When taking this limit, one has to exercise some care in treating 
the singular functions and distributions in Eq. ( [4.27] ). It follows from the identity ( |4.30|) 
that for any test function f{x) we have 

T{x; e) f{x) dx = T^'\x] e)Q{l - x - tx) [f{x) - /(I)] dx 
Jo Jo 

= I T^''\x] t) \f{x) - /(I)] dx + /'(I) / T'^'\x] e){l-x)dx + .. 

JO "'l/(l+<:) 

(4.36) 

Defining '+' and '++' prescriptions such that 

F{x)+ = F{x) - 5{l - x) F{z) dz , 

. (4.37) 
F{x)++ = F{x)+ + 6'{1 - x) / F{z){l -z)dz, 

JO 

and recalling that 

f'5'il-x)fix)dx = f'il), (4.38) 

JO 



OK 



we should therefore define the singular terms in Eq. ( |4.27| ) at small e, up to terms of order 
e, as follows: 

^ ^ +e6'{l-x) 



1 — X {1 — x) 
ln(l — x) I ln(l — x) 



1 — X \ \ — X 

e e 



{l-xf {l-x)\^ 
-e5'{l~x) 



+ (elne-e) 5'{l-x) 

+ 

elne(5'(l — x) 



(4.39) 



The small-e behaviour of JF is thus of the form 

J^[x] e) = -P{x) In e + c(x) + ^(x) e In e + e + O (e^ In e) , (4.40) 

where 

Pi.) = _ (1 + .) ^ ,(1 _ = (i±^) ^ (4.41a) 

,M . 2( '"l'\-f-'l )^-(l + .)ln(l^j + l + 4. 

--(5(l-x) (4.41b) 



2(l-a;)+ 4 
4 

c,(a;) = -^^ — — 2(2 + X + 6x2)+ 9 5(1 -x)+5'(l-x) (4.41c) 

Mx) = +2(2 + x + 6x2)lnfl^ 

\1 — x/^ Vx"' 

^ + 7^ — + 8 + 13x - 16x2 - 8 5(1 - x) - 4 5' {I - x()4.41d) 



The corresponding formula in moment space is 

J^iv(e) = -Pn^'^e + cn + Qn e\ne + hN e + O (e^Xne) , (4.42) 



where 



= + (4.43a) 

31 1\^ ^2-^92 3 1 1 



4 2 12 

gj^ = N + 8- — - — — ASi (4.43c) 

/ 4 2 1 2 \ 

= (^ + «-iv-ivTT-ivT2)^-^^? + «^^ 

-5^-3 + ]^ + ]^ + ]^+(FTT)J+(]^ (4-43d) 



with 



^-1 1 „2 

'^.= E-^; 5i = ^(iV)+7£; = lniV + 0(l/iV) , ^2 = - - ^'(iV) . (4.44) 
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Figure 6: Coefficient of higher-twist contribution to DIS. 



Power corrections. From ( [4 .401) and (|4.42| ) we can now evaluate the dispersively- 
generated power corrections to the usual (logarithmic) scaling violation using Eq. ( [4.24| ). 
The leading non-analytic term of T{e^ N) for small e is g^^ In e, leading to a power-behaved 
contribution of the form 



d 



rflnQ2 



In Fr{N,Q') = gn^^ 



(4.45) 



We shall not undertake detailed phenomenological studies in this paper, but only indicate 
how the above result could be applied. In the analysis of DIS data, the non-perturbative 
contribution is normally parametrized as a simple multiplicative correction of the form 



1 



C{x) 



(4.46) 



The 'higher-twist' coefficient C{x) is then found to be a steeply increasing function of x, 
as illustrated by the data points in Fig. ^ ^ . 

In our approach, the leading non-perturbative contribution is of the form 



Fr{x,Q') = -A',g{x)®Fr{x,Q')/Q' 



(4.47) 



where '®' represents the convolution corresponding to a product in moment space. Thus 
we predict a coefficient in Eq. ( [4.46| ) with a weak (logarithmic) dependence. 



-A',gix)^Fr{x,Q')/F^ 



(4.48) 



The curve in Fig. ^ shows the prediction of Eq. ( |4.48| ) at = 10 GeV^, assuming 



A2 = —1 GeV . For F2 we have used the valence part of the MRS A parametrization 



3M . We see that the form of the observed higher-twist correction is well reproduced. The 



07 



steep rise at large x results mainly from the most singular term 6'{1 — x) in g{x), which 
generates the behaviour 



C(x,Q2)~A;^lnFr(x,Q^) (4.49) 

at X 1 . However, the less singular terms also play a significant role in the region of x 
shown in Fig. ||. 

4.3 DIS sum rules 



Adler sum rule. As we have already mentioned above, due to the identity ( [4.30[) the 
Adler sum rule (for a quark as a target) 

dx F2{Q^ , x) = I (4.50) 
does not acquire either perturbative or non-perturbative correction to this order. 

Gross-Llewellyn-Smith sum rule. This sum rule concerns the first moment of the 
F3 structure function. For a target hadron h one has 

G{Q^)= [' dx[F^\Q^x)+F^\Q^,x)] = Go dxFsiQ'',x) , (4.51) 
Jo Jo 

where Go = 6-B — 25", with B and S the baryon number and strangeness of the hadron h, 
is the parton-model value. Here F3(Q^,x) stands for the non-singlet structure function 
for lepton scattering by a single quark. 

According to our procedure we first compute the matrix element squared for massive 
gluon emission corresponding to F3. A simple calculation gives (apart from the coupling 
as and the colour factor Cp/^n) 

Md = M3-M2 = — ^2 (4-52) 

{y - xe) 

where A^2 is given by ( |4.25| ). Since the virtual contributions to F2 and F^ are identical, 
the difference between their characteristic functions is given by 

Tsix, e) - J^2{x, e) = /' dy Md{x, y, e) = -2 (9 xe - e - 4) xh In ( ^ ) 

Jxe/{l-x) \[1 — Xe)[l — X) J 

- (1 + x) h e 18x^e+ 18x^ + 7xe - lOx + 5e - 4 ^ H . 

1-xe V (1-x)^ 1-x y 

(4.53) 

This difference, determined by real emission only, is defined for e < (1 — x)/x and vanishes 
at the phase space boundary as 

J^3{x, e) - J^2{x, e) = -{l + x)[ + . . . . 

\1 — X X/ 

The first moment reads 

nl/(l+e) 



J^3(e, A = l) = / dxJ^s{x) 

4^)-^-(^4)(-^Ni..)-.'H-i). 

(4.54) 
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It has the following small-e behaviour, relevant for power corrections: 

83 



^3(e,iV = l^ 
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(4.55) 



Thus the one-loop correction to the first moment of the quark structure function F3 is 



Fs{Q',N=l) = -^ 



27r Jo /i^ 



«efr(/^')^3(e,A^ = l), 



where 



T,U.N = 1] 



2 3 V3 3/ e 3e^ 



(4.56) 



(4.57) 



As usual (for a collinear-safe quantity) this function vanishes for both large and small e. 
The area under the curve gives the well-known one-loop perturbative contribution 



Gn 



3Cf 

4:71 



-a 



(4.58) 



We may now evaluate the dispersively-generated power corrections, which are determined 
by the small-e behaviour 

j='3(N = l,e) = e f-lne+ — ) +e^ f-lne- — ) + regular terms. (4.59) 
V3 9 / V3 36/ 

The two non-analytic terms generate power-behaved contributions of the form 



° l3g2 ^ 3g4 



(4.60) 



The estimated coefficient of 1/g 13^ used in a recent analysis at g ~ 3 GeV by the 
CCFR collaboration ||36| is —0.09 ±0.05 GeV^. This implies a small negative value of A2. 



However, at such a low value of the l/Q"^ contribution could also be significant. 
4.4 e+e~ fragmentation function 

Our procedure for e"^e~ annihilation is similar to that presented for deep inelastic scat- 
tering. The e"'"e~ fragmentation function F{x, Q^) is given in terms of the characteristic 
function 

T{x; e) = T^''\x; e)Q{l - e - x) + H(e)(5(l - x) (4.61) 
where J-'^^'^ is the real contribution from annihilation into quark-antiquark-gluon. 



X, e 



1 — X 



(1 +x) -2e 



In 



(x + e)(l — x) 



1 3 + 4e + 3e^ e(l + e) e^ 



(4.62) 



+ 



X + e 



and Vt{e) is the timelike virtual correction given in Eq. (|4.3b|) . In this case we do not 
have an exact sum rule of the form (^4.30| ): instead 



Vi(e) = 1-1' T^'\x- e)e(l - t - x) dx + O ie' 
2 Jo ^ 



(4.63) 



on 



However, this is sufficient for our purposes, because it means that for any test function 
f{x) we can still write 

/V(x;e)/(a:)rfx = ^/(l) 
JO 2 

+ f :F^^\x-,e)[f{x)-f{l)]dx + f\l) C ^M(x;e)(l-x)dx + ... 

(4.64) 

up to terms of order e. Using the dictionary of singular terms ( [4.39| ), we thus obtain in 
this case the small-e behaviour ([4.40|) with the same splitting function P{x) but with c(x), 
g{x) and h{x) replaced by 

= 2(^i^^^^^^^ _(l+a;)ln[a;(l-x)] + i(l+x) 



2(1 4 

~g{x) = - ^ +2 + 5'{l-x) 
U — x}^ 

~ , , /ln[a;(l -a;)l\ , , , 1 2 

h{x) = ^ _21nxl-x +7^ ^ + - 

\ 1 — X ) ^ (1— X 

-h5{l-x)-A5\l-x). (4.65) 
The term 2/x in /i(x) needs some interpretation. It results from setting 

In fl + + ^— ~ - (4.66) 

V xj X + e X 

which is not valid unless x S> e. Hence the expression for h{x) is not to be used in the 
small-x region, or for moments that are sensitive to it. In particular, in moment space we 
have a result of the form (|4.42|) with J^at = | + C (e^) for iV = 1, while for > 1 

31 1\^^9^93 1 2 a 



~9n = N-l + ^ + ASi (4.67b) 



hN = [^N-l + -)S^ + 2St + 6S2-5N + j;^ + — . (4.67c) 

For phenomenological applications, taking the same approach as for DIS, we may 
parametrize the leading power correction to the fragmentation function as 

Fix, Q') = F^\x, Q') {l + (4.68) 

where 

C{x, Q2) = ~g{x) ® FPT(x, Q')/FPT(x, Q^) . (4.69) 

Fig. 1^ shows the resulting prediction using the same parameter value as for DIS, 
= —1 GeV^. For the perturbative contribution F^^ we used the parametrization of 
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Figure 7: Coefficient of power correction to e^e fragmentation function. 



the light quark fragmentation function at Q = 22 GeV given in Ref. ||3^. The coefficient 
function C is larger than that for DIS, but qualitatively similar in form. The values of 
in e^e~ annihilation being much larger, the predicted power correction is probably too 
small to be detectable. 

In Ref. the non-perturbative contribution to the fragmentation function was 

parametrized small shift in the value of x, 



X 



X 



(4.70) 



where = —0.14 ± 0.10 GeV. However, a parametrization with in the place of 

1/Q was also found to be acceptable, the magnitude of the correction being practically 
consistent with zero. A small shift in x is clearly equivalent to a correction proportional to 
F', the x-derivative of F, as expected from the 6' term in g{x). As shown by the dashed 
curve in Fig. 0, a correction of this form is similar to our prediction at large x. The 
dashed curve corresponds to C{x) = —0.35F'/F, which, at 22 GeV, would be equivalent 
to ho ~ -0.016 in Eq. ( ^70[) . 



4.5 Event shape variables 

The analysis of power corrections to event shapes in e^e~ final states proceeds in the same 
way as for hard cross sections, except that in this case the observables are constructed in 
such a way that there is no virtual contribution. Thus for the mean value of some generic 
shape variable y, defined so as to vanish in the two-jet limit, we have the characteristic 
function 

= jrM(g) = J dx dx Meeix, x, e) y{x, x, e) e[{l-x){l-x)-e]e[x+x-l+e\ , (4.71) 
where Aiee is the matrix element given in Eq. ( [4.6| ). 



Q1 




Figure 8: Derivative of characteristic function for mean value of thrust. 



Thrust. In the case of the thrust variable T, for example, we define y = 1 — T, which 
vanishes in the two-jet limit, and 



y{x, X, e) = min{(l - x), {1 - x), {1 - y'x^ - 4e)} . (4.72) 
where As before, the quantity to be inserted in the fundamental equation 



( p.5|) is JF, the logarithmic derivative of the function ([4.71|) . 

We see from Fig. |^ that JF for the mean thrust decreases much more slowly at small e 
than the corresponding function for the quantities studied earlier. In fact the behaviour 
at small e is 

T{e) ^Ayfe + 0(eWe) , (4.73) 

which is shown by the dashed curve. Thus the leading term is non-analytic and of square- 
root type. It follows that the leading power correction to (1 — T) will be of order as 
observed 

It is easy to see how the y/e behaviour arises. It comes from the contribution to 
the derivative from the first theta-function in Eq. ( ^.711 ), i.e., from the soft phase-space 
boundary, which gives 



jF(e) ~ ey (ixfix A1ee(2^, 2:, 0) min{(l — x), (1 — x)} (5[(1 — x)(l — x) — e] 6[x + X — 1] 
~ 4 dx = 4v^ . (4.74) 



Notice that we obtain in this approximation 



(1 — T)~ — / dx dx aeff[{l — x){l — x)Q'^] Meeix, x,0) mm{{l — x), {1 — x)} , 
2n Jo Ji-x 

(4.75) 
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which corresponds to the result obtained in Ref. 0: the argument of aeff is the maximum 
gluon virtuahty, which is equal to its transverse momentum. 

In terms of the non-perturbative parameters ( [4.1D , the power correction to (1 — T) is 
given by the V = \ moment of Sa^a as 

(1_T)^P^4^. (4.76) 

The analysis performed in Ref. showed that this gives a good description of the data 
for 

/li~0.25GeV. (4.77) 

In Ref. 0], the same 1/Q dependence as in ( |4.76| ) was obtained, but with a different 
coefficient. The difference arises from the normalization factor in the definition of the 
thrust: we have normalized to the sum of the final-state energies, whereas in the sum 
of momenta is used, corresponding to inserting a factor of 



2/ (2 



Xg + Jxl-At) (4.78) 



into Eq. (f4.72|) . Even though the thrust as usually defined is normalized to the sum of 



final-state momenta, the factor (|4.78|) should not be included, because it corresponds to 



finding a real massive gluon in the final state. The massive gluon in the calculation always 
decays into massless quarks or gluons, and so the sum of the final-state momenta should 
be set equal to the sum of the energies. 

As mentioned in Sect. 3, the dispersive method applies directly to quantities that are 
fully inclusive with respect to gluon branching. Then the entire effect of branching is 
to make the coupling run. However, event shapes are sensitive to the structure of the 
final state, and branching may lead to a different value of the observable. In the case of 
thrust, for example, the value is unchanged only if the products of gluon branching fall 
into the same hemisphere, which includes quasi-collinear branching and therefore gives 
the dominant contribution to (1 — T) ~ «cfr/7r. Branching into opposite hemispheres, 
corresponding to a genuine four-parton contribution to the thrust, gives a correction 
~ (ttefr/Tr)^. Both give rise to 1/Q power terms |2^. The extent to which terms of 



higher order in cteff affect the magnitude of the power-behaved contribution remains to be 
established phenomenologically. The result ( [4.77| ) suggests a typical value of the effective 



coupling in the small momentum region such that acfr/7r ~ 0.2. This gives grounds for 
optimism that higher-order terms may be controllable. 

Similar results to ( |4.76| ), with different coefficients of Ai/Q, may be obtained for the 



mean values of a variety of e^e~ and DIS final-state event shapes P,^. 

Since the presence of a 1/Q correction introduces a large non-perturbative contribu- 
tion, one would prefer for some purposes to define event shape variables for which the 
predicted coefficient of 1/Q is zero. Such variables should be more suitable for testing 
perturbative predictions and for measuring a^. As suggested in one can find linear 
combinations of variables such that the predicted 1/Q terms cancel, at least in the mean 
value. In addition, there are shape variables for which such terms vanish because the 
small-e behaviour of the characteristic function is not of the square-root type. 

Three-jet resolution. An example of a shape variable without a leading-order 1/Q 
correction is the mean value of the three-jet resolution variable 1/3, defined according to 
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Figure 9: Derivative of characteristic function for mean value of y^. 



the Durham or k± algorithm In lowest order, in the region x > x, we have 

?/3 = mm{yqq,yqg} 

where 

ygq = (1 - cos 6*55) 

1 r;^.2 ^2 



yqg = 2 min{x - cos6'qg) 



with 



cos 6', 



99 



COS 6, 



19 



1 - - e 

XX 

Xg — 2(1 — X — e)/x 



x^ - 4e 



(4.79) 



(4.80) 



(4.81) 



The resulting behaviour of JF for the mean value of y^ is shown in Fig. ^. The small-e 
behaviour in this case is 

J^(e) e\n\ + 0{e\ne) . (4.82) 

It follows that the leading power correction to the mean value of ys should be of order 
In Q/Q"^. We find numerically that the non-leading logarithms are such that 

J^(e) ~ e(ln2e + 31ne + 4) , (4.83) 

as shown by the dashed curve. In terms of the moment integrals ( [4.1|) , we therefore expect 

{ysf' ^ -i2A', InQ' - 3A', - A'^/Q' . (4.84) 
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If the three-jet resolution is defined instead according to the JADE algorithm, then 
Eqs. (|4.80| ) become 



yqq = ^xx{l - cos9qg) 

Uqg = Xg {1 - cos 6qg) . (4.85) 

In this case the behaviour at small e is of square- root type, leading to a 1/Q power 
correction. This is probably why the Durham algorithm has been found to require smaller 
non-perturbative corrections . 

A final point to be noted from Figs. |^ and ^ is that the characteristic scale of hardness 
for event shapes, as typified by the peak of the characteristic function, lies far below Q^, 
at ~ 0.05Q^. This is a general feature of quantities which probe final-state structure 
in time-like processes, such as event shapes and e+e~ fragmentation functions, and could 
explain why fixed-order fits to such quantities favour small scales and large values of the 
coupling. 

4.6 Drell-Yan process and K factor 

The matrix element squared for the Drell-Yan (DY) process with emission of an off-shell 
gluon (yU^ > 0) is given by 

g + r^g + 2 re - 2 r^a^ + 2 + eV^a 

Mdy = 2 "2 , (4.86) 

r {e + a) 

where e = /x^/Q^, r = Q'^/s and a = q'^/Q'^ with the momentum integration variable. 
The phase space integration is given by 

= da-^ e{R) , 
v R 



with 

R = A-4:T^a, A = l-2T + r^ -2Te-2r^e + eV^ . 
The real part of the DY characteristic function is then 



(e, r) = rf$ Mdy = -4 + 4 tanh-^ ( ^ ^ \ 

JO \ \ — T — Te ) 1 



r — re 



(4.87) 



We introduce the "rapidity" variables 



c=l±i (4^ 
2(coshr/ + c) ' 2V^ ' ^ 



with the phase space 



c -|- \/(? — \ sinh "qdr] c -|- \/ — \ 

dr = — _ _ 



2 (cosh r^ + c)'' "''^ 2(1 + c) (1 + V^)^ ' 
One then has the simple form 

cosh 7] + c \ cosh r] J \ (cosh t] + c) J 



The complete characteristic function is obtained by including the virtual contribution. In 
moment space one has 

r-M r /"'"max 71 r I \ 

(e) = / dr r^-i ^^i^(e, r) + H(e) . (4.90) 

This function is given to order e in Appendix A, Eq. ( [A.6|) . The DY distribution moments 
are then given by 

<m = ^rK «eff(/^^) (e) , (4.91) 

/vr JO /i 

and the moments of the Drell-Yan K factor are 

^iv(Q') = ^r% «eff(/i') t^{e) , (4.92) 

ZTT JO yU 

where the characteristic function is 

/C^(e) = J-^^(e)-2^^^^(e) (4.93) 

with jF^^'^(e) the corresponding moment of the DIS characteristic function ( [4.26[ ). As is 
well known, the combination ( |4.93| ) is coUinear safe since the collinear singularities in the 
Drell-Yan and DIS terms cancel. Its logarithmic derivative K,n{^) vanishes both for large 
and small e, as shown for the first four moments in Fig. ^ At small e, K,n{^) has the 
form 

/Cw(e) =CN + e{fN \n\ + gN\ne + hN) + (e' In' e) , (4.94) 

(see Appendix A, Eq. ( |A.7| )), which shows that, apart from logarithmic enhancement, the 
power correction to the Drell-Yan cross section is of order as in DIS, in agreement 

with the result of 0. 

The expansion of 1Cn{^) at small e has the form 

/Cjv(e) = e [/^(Ine + 2Sif + ^jv(lne + 2^i) + h^] {l + O (e)) , (4.95) 
where at large iV, corresponding to r ^ 1, we have 

/n-- N , 2iV2 , fiN ~ AN"^ . (4.96) 

Thus the expansion parameter at large becomes eA^^, and for /Cjv one has approximately 

M.)^2eA'^(.n^ + 2)(l + o(!!i(lp)). (4.97) 

where we have used 

lne + 2Si ~ ln(eA^2) . (4.98) 

Correspondingly, at large r the power correction to the K factor may be expressed in 
terms of the non-perturbative integrals ( [4.1|) as 

This means that the coefficient of the power term increases much more rapidly as r — > 1 
than does the DIS coefficient C{x) as x — > 1 in Eq. ( [4.46[ ). 




Process 



Quantity 



Power 



Coefficient (ffist order) 



Eq. 



e ' e 
r 
DIS 
DIS 



e+e 



DY 



R 

F2 
GLS 
F 

(1-T) 
K 



-UA'qL - 12A'q + 6A'^ 



4 A' 
3^2 



-2A'2L + 3A'^ + A'^ 
2fA',L - (2/ + - M' 



(EH) 
(14:211) 

O) 
HI) 
(PI) 

iH) 

(EE) 
(B 



Table 1: Summary of power-behaved contributions. 
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4.7 Summary 



We summarize in Table 1 the results obtained for the leading power-behaved corrections 
to the quantities discussed above. In each case the predictions are expressed in terms of 
the effective coupling moment parameters given in Eq. ( |4.1| ), and L = lnQ^ (Inm^ in the 
case of Rr). 

At present no systematic phenomenological tests of these predictions have been per- 
formed and the values of the moment parameters are correspondingly uncertain. The 
parameter Ai specifies the lowest-order 1/Q contributions to a wide range of event shapes 
in e~^e~ (and DIS) final states, of which we show thrust only as an example. A value of 
Ai ~ 0.25 GeV is suggested by the 1 GeV/Q power-behaved contribution seen in the data 
on (1 -T). 

As discussed in Sect. 4.2, the l/Q^ contributions to the structure function F2 in 
DIS appear consistent with a log- moment parameter value ~ — 1 GeV^. The power- 
behaved contribution to the Gross-Llewellyn-Smith sum rule for F3 may suggest a smaller 
value. In either case the corresponding contribution to e+e~ fragmentation functions is 
probably too small to be observed. The other 1/Q^ effects listed, for the mean three- 
jet resolution in e~^e~ final states and the i^-factor in the Drell-Yan process, show a 
logarithmic enhancement and therefore involve also the additional log-moment parameter 



5 Conclusions 

In this paper we have started to explore the possibility that the concept of the QCD 
running coupling can be extended, in a process-independent way, down to small momen- 
tum scales, at least in an effective sense. At low scales confinement physics dominates, 
undermining the very possibility of applying the language of quark and gluons, the only 
language we are able to use at the present level of our understanding of QCD. The hope 
is that, in spite of all the richness and complexity of the spectrum of hadrons, the bulk 
of the non-perturbative physics will reveal itself in a smooth way in sufficiently inclusive 
observables, and that it can be taken into account by extending the notion of the pertur- 
bative QCD coupling. The coupling aes thus defined is intended to measure the effective 
interaction strength even at large distances, allowing the extension of the quark-gluon 
language beyond its original domain of applicability. 

In an inclusive observable F, confinement reveals itself via departure of the measured 
value from the perturbative prediction F^"^ at decreasing values of the relevant hard scale 
of the process, Q^. At sufficiently large Q^, F^"^ varies slowly (logarithmically) with 
Q^, since it is given by a finite-order expansion in ^^(Q^). The departure is typically 
power-behaved (see Table 1) 

F^P = F(g2) - F^'^m ^ [C^A,, + C^A',^ + CsAl] ■ Q-'^ , (5.1) 

where the dimensionless coefficients Ck are at most linear in InQ^. In this paper we have 
developed a method of calculating the pro cess- dependent exponents p and coefficients Ck 
for a wide variety of hard observables which do not involve gluons in the Born approxi- 
mation (quark-dominated processes). The dimensionful parameters A2p, and are 
universal and are given by (log-)moment integrals of the effective coupling modification 
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5aes in the small momentum region. These are new phenomenological parameters which 
should be studied experimentally. We have computed the corresponding power-behaved 
contributions to first order in da^s . 

Our method is based on using a dispersion relation for the running coupling a^- The 
dispersive machinery then leads to the basic representation ( |3.5| ) for the observable F in 
terms of an integral over /x^ of the product of a characteristic function and the effective 
coupling Q;efr(/i^). The latter is related to the standard as by the operator equation 

2n sin(7rP) ^ 2d 



This means that ttefr and are practically equivalent in the perturbative region. At low 
scales, aeflf and as may differ substantially from each other and from expectations based 
on the perturbative /^-function. However, we have argued that the non-perturbative 
modifications to aesifi"^) at low /i^ should be such that no power corrections to are 
generated at higher scales. It follows that the /z^^-moments of these deviations have to 
vanish for the first few integer values of p. Thus in order to generate power corrections, 
JF must contain terms which are non-analytic in /i^. 

The characteristic function JF is specific to a given observable and can be analysed 
by Feynman diagram techniques. We have computed the characteristic functions for a 
number of observables in various processes to first order in tteff- In this order, ^ is a 
function of the ratio e = /i^/Q^ and is given by the one-loop diagrams with a gluon of 
mass equal to the dispersive variable /x, < /i < 00. The leading non-analytic term in 
the e — s> behaviour of determines the power p and coefficients Ck in Eq. ( [5.11 ). 

Let us stress again that we work with a massless gluonic quantum field. However, in 
the characteristic function the dispersive variable /i enters as a gluon mass in Feynman 
denominators and the phase space. The role of a small gluon mass as a trigger for long- 
distance contributions to hard processes has been recognized and exploited in the recent 



literature, see [Pl,f7,12,17 



An attractive feature of the dispersive method is that for a given process it suffices 
to compute a single function to obtain power corrections to all associated observables. 
For example, the matrix element squared for e~^e~ annihilation into massless qq and a 
massive gluon gives the dispersively-generated power corrections to Re+e-, to the non- 
singlet fragmentation function, and to all event shapes (thrust, y^, etc). The results 
obtained for these and other hard process observables are summarized in Table 1. 

As emphasised in Table 1, our calculation of the coefficients of power-behaved terms 
is first-order in the effective coupling a^s- The accuracy of the first-order estimate will be 
reasonable if aeff/n happens to be numerically sufficiently small in the important kinematic 
region. This is a question that can only be answered by comparison with experiment. The 
same qualification applies even more strongly to event shapes which, by construction, are 
sensitive to final state structure and not fully inclusive. Multiparton contributions to such 
quantities are higher-order in a^s, and may or may not be suppressed, depending on the 
typical value of aefr/vr. 

On the topic of higher-order contributions, we would like to end with the following 
remark. A characteristic feature of the leading-order power-behaved terms that we have 
calculated is that they are enhanced near the boundary of phase space, where the addi- 
tional gluon is soft. One thus observes the following general pattern of enhancement of 
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power terms and the corresponding phase space boundaries: 

DIS coefficient functions = ^ (1 - x)Q^ ~ 

Q 

DY ir-factor = — A^ ^{l-rfQ^^K^ 

Q 

differential thrust distribution =- A^ =^ (1 — T)Q'^ ~ A^ 

(1 - T)Q^ 

etc. 

The power contribution becomes of order 1 at the edge of phase space, where the invariant 
mass of the final-state hadronic system is squeezed down to a finite value ~ A^, so that 
the process becomes quasi-elastic and can no longer be treated as hard. As has already 
been remarked 0, it is possible that such effects, being closely related to the Sudakov 
form factor, can be shown to exponentiate and factorize universally. 
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Appendix A : Drell-Yan cross section 

The real contributions of the moments of the Drell-Yan characteristic function are given 
by 



[ 1] sinh rj 



2c 



N 



+ 



drj 
7] sinh 7] 



[ coshr/ (cosh ?7 + c)^ coshr/ (coshr/ + c)^+2 
Introducing 



— sinh^ r] 



(cosh T] + c) 



N+2 



hie) 



drj 1] sinh rj 



In^ (c + V^^^ 



hie) 



cosh?7(coshr7 + c) c 

dr] sinh^ r] _ 2 In (c + y/c^ - 1 
(coshr^ + c)^ (c^ — 1)^/2 



TT 

T 



2c 



C2 - 1 ' 



we obtain 



'cW?3I' 
2c 



N 



_Dn^i + Dn+i]Ii{c) + 



A^(A^+1) 



Dn-iHc) 



(A.l) 

(A.2a) 
(A.2b) 

(A.3) 



with 



-"+1 / d 



n\ \ dc J 

For the purpose of extracting the first power correction, the expressions ( A.2 ) may be 
expanded as 



hie) 



hie) 



21n^(2c) ^ TT^ _ ln(2c) ^ ^ , 
c 2c c'^ ^ 



2^2 ln(2c) 2 



c 



(A.4a) 
(A.4b) 



The operator Dn acts as follows: 



r(fc + ^) k+i 

T{k)n\ 



D„ { ln(2c)c-^' } = Ei^±^c-'=+i(ln(2c)+^(A;)-^(fc + n)) 



Dn{\n\2c)c-''} 



T{k)n\ 
T{k + n) / 



T{k)n\ 

Applying the above results, one arrives at, to first order in e. 



c-''+^ {ij'{k + n)- i:\k) + [ln(2c) + ^(fc) - ij{k + n) . 



M 



(0) 
N 



(Ine + 2SiY + ^^ + 2(lne + 2^i) 



1 1 

TT + 



452 



(A.5) 



M^^^ = -Ar[(lne + 2Sif + n'] + 2(lne + 2^0 



+ 3A^ - 3 + 



1 



6A^(A^ + 1) + 4:NS2 + 16 



1 



A 1 



Taking into account the virtual correction, one obtains 



+ e |-(iV + 2)(lne + 2S^f + 2(lne + 2Sx 
-6N(N + 1) + 4 (NS2 - 2Sl) - (at - ^) + 12 - 



(A.6) 
Constructing 

IC^ie) = ^^^(6) - 2^^^^(6) , 

one derives the following expression for the characteristic function of the i^-factor in 
moment space: 

/Cjv(e) = cn + e{fN In^ e + ^jv In e + /ijv) 

where 

^,^2 r.x/o2 2\^47r2^4 6 2 
c. = 2(^? + ^2) + (-3+- + ^)5, + - + l---^ + ^ + 

/jv = -(iV + 2) 

^AT = 2Ar(Ar + 2) - 4(A^ - 2)Si -22 + 2 



(TV +1)2 

19A^2 ^ 30iv + 8 
N{N + 1){N + 2) 



- eiV^ + 4Ar + 18 - AT - - TT 



= -4(7V - 1)5,2 + (47V2 + 107V - 28 + 8 ^^j^' j 5^ + 4(7V - 3)^2 

9N^ + 61N^ + 122iV3 + 104A^2 ^ 43^^ ^ 

^ A^2(Ar+ l)2(iV + 2)2 ■ 

(A.7) 

Its logarithmic derivative may be written in the following form, which is suitable for 
studying the large A'" regime: 

JCNie) = e [/^(In e + 2Sif + g^iln e + 25i) + hj, 

where 

/„ = JV + 2 



Ar(Ar + l)(A^ + 2) y 

/ 4 2 12 \ 

hN^-2 (^ + 16-^-^-^) 5i + 4(Ar - 1)2 _ 4(Ar - 3)^2 + 20^2 

, 2\ 2 ^5N^ + 1SN'^ + 7N^ -lON"^ -16N -S 



3 J N^{N+iy{N + 2)^ 

(A.8) 
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